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Abstract. We investigate the evolution of the interaction of the magnetoplasmon 
resonance with the harmonics of the cyclotron resonance as the confinement of an 
electron gas in a quantum wire increasingly deviates from the parabolic case. The 
occurrence of the Bernstein modes is observed in a time-dependent Hartree model of 
a two-dimensional electron gas in a single quantum wire. 



1. Introduction 

Quantum wires are interesting for probing the electron-electron interaction in the 
transition regime from quasi two-dimensional systems to quasi one-dimensional systems. 
These systems can be fabricated starting from two dimensional semiconductor 
microstructures (2DES), for example in AlGaAs-GaAs heterostructures. The wires may 
be created by a gate structure or by etching techniques. The presence of a magnetic 
field is manifest by many interesting quantum features, most notably the quantum Hall 
effect and the fractional quantum Hall effect. Far infrared response is a powerful method 
for determining the charge-density spectra in the transverse direction. For transverse 
parabolic potentials, the motion of the electrons may be separated into center-of-mass 
motion and relative motion. The center-of-mass motion is decoupled from the relative 
motion of the electrons, and takes place in a transverse potential with the effective 
frequency (lUq + uj 2 ) 1 ^ 2 , where Tiujq is the confinement energy and to c is the cyclotron 
frequency. The time- dependent external potential <p e representing the long wave length 
far-infrared radiation only couples to the center-of-mass motion causing a single peak 
in the FIR spectrum at the plasmon frequency VL P = (cu 2 + uJq) 1 ^ 2 equal to the effective 
transverse frequency |T], Therefore, according to Kohns theorem^], although a rich 
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spectrum would be expected, there is only a single resonance. For small systems with 
general confinement potentials deviating strongly from the parabolic case, the single- 
particle excitations may overlap with the collective excitations causing a fine-structure 
in the plasmon peak[]3], |J. Any deviation from parabolic confinement will couple the 
internal and collective motion. 

Bernstein showed that in 3D homogeneous systems the dispersion of the 
magnetoplasmon as a function of the magnetic field will have anticrossings at finite wave 
vectors (k _L B) around Q p = nu c where n = 2, 3, • • •. The calculation was based on 
solving the classical Boltzmann equation and the Maxwell's equations self-consistently. 
The strength of the anticrossings depends on temperature, wave length and magnetic 
field in a complex way and decreases for the higher harmonics. An interaction between 
the magnetoplasmon resonance and the harmonics of the cyclotron resonance, a so-called 
Bernstein modeg, have been observed in 3DES[|, 2DES0 ||, IDES (wires) and ODES 
(dots) [|J . For confined electron systems in reduced dimensions the breakdown of Kohns 
theorem in nonparabolic potentials is sufficient to observe the Bernstein modes instead 
of considering absorption at a finite wave vector. 

We study the far infrared absorption in the integer quantum Hall regime. Here 
we consider a broad range of deviations from parabolic wires in the time-dependent 
Hartree approximation. Increasing deviations give richer excitation spectra when the 
magnetoplasmon resonance starts to couple to increasingly higher harmonics of the 
cyclotron resonance. We systematically investigate how the splitting and its location 
depend on the strength of the deviation from the parabolic confinement. 

2. Model 

We consider a strictly two-dimensional electron system lying in the x-y plane. The 
motion in the z direction is neglected since the electrons are confined to the lowest 
subband at the low temperature attained in experiments. We use the Hartree 
approximation (HA) to reduce the many-particle Hamiltonian to a single-particle 
Hamiltonian for each electron in an effective potential approximating the electron- 
electron interaction. For the FIR absorption we use the corresponding time-dependent 
approximation describing the self-consistent linear response of the 2DES to an external 
homogeneous time- varying electrical field (the random phase approximation, RPA). In 
a constant perpendicular magnetic field B the cyclotron frequency and the magnetic 
length are uj c = eBj (m*c) and l c = [%/ (m*c<; c )] 1//2 , respectively, where m* is the effective 
mass. The dielectric constant of the surrounding medium is noted by n. It is convenient 
to introduce the constant perpendicular magnetic field with the vector potential in the 
Landau gauge A{f) = (—By,0). 

The effective Hamilton operator for a single electron in the confining potential V c (y) 
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is 

h - ( v2 - 14) + \ m '^ + VM + ^ • « 

where V#(r) is the self-consistent Hartree potential representing the direct Coulomb 
interaction between one electron and the total charge density of the 2DES. Periodic 
boundary condition in the longitudinal direction of the wire gives a Bloch-type single- 
particle wave function 

q nk (x, y) = -=e ikx ip nk {y) , (2) 



J x 

where the longitudinal wave vector k = p ■ 2tt/L x with p G Z and n is the transverse 
quantum number. The length of the wire, L x , is assumed to be much larger than the 
effective width of the wire, therefore the Hartree potential only depends on the transverse 
coordinate [H 
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with L given below. It has been assumed that a neutralising background charge 
exists. For noninteracting electrons in a simple one-dimensional parabolic potential 
V(y) = |m*CG>Q?/ 2 , the single-particle eigenf unctions are 



4>nk(y) = 7= , — l -^ H n 
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where y k = k ■ + It) is the center coordinate, which generally does not equal 

the expectation value of the transverse coordinate j/flKj. The nth Hermite polynomial 
is denoted by H n . The confinement length is defined as l = \h/ (rrfujo)] 1 / 2 . The 
electron is localised within the effective length L = [/i/(m*f2)] 1//2 , replacing the magnetic 
length l c , and the cyclotron frequency uj c is replaced by an effective frequency defined 
by Vt = (lUq + cu 2 ) 1 / 2 . The eigenenergies corresponding to the eigenstates (§) are 
E nk = hQ(n + 1/2) + (h 2 k 2 /2m*) x + In the case of vanishing confinement, 

(2DES) Iq — > oo, the energy bands reduce to the familiar dispersionless Landau levels. 
The effective single-particle Hamiltonian is diagonalised using the wave functions of the 
noninteracting electrons (£|) as a functional basis and the self-consistent solutions are 
obtained by iteration. 

The response to a time-dependent perturbation may be found by the random phase 
approximation where exchange and correlation effects are neglected. In this mean- 
field approximation the noninteracting Hartree quasi-particles move in a self-consistent 
potential given by the external perturbation and the response of the charge density of 
the electrons n s . The selfconsistent potential, (j) s (f,u), may be found from the external 
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potential 4> e (r, to) 

(a\M«>m = (<*\M<*)\P) + E ^^4 (^)17) , (5) 

^ ft^ + (e 7 - e S ) + i0 + 

where a denotes the longitudinal and transverse quantum numbers n and k. f a is the 
Fermi occupation factor and the Hartree matrix elements are defined in terms of the 
Hartree ground state wave functions 

H a ^ 5 = -JdrJdr ^-^ . (6) 

A method to find <j) s based on a Fourier transformation of the position coordinates and 
repeatedly solving (|IJ) for all the needed values of uj is shown in Ref. [ffl. Here we 
present an alternative method transforming ([5]) such that (j) s {oS) can be simultaneously 
obtained for all values of uo from a matrix eigenvalue problem. The latter method 
considerably speeds up numerical calculations and we have checked that they deliver 
exactly the same results. The equation for the selfconsistent potential (HD simplifies 
for transverse external fields. The longitudinal quantum number k is conserved in the 
transition, so that each single electron transition may be labeled by k and two transverse 
quantum numbers n and m. By applying a time dependent but spatially constant electric 
field, E e (f,t) = —yS e exp(—iujt), the Fourier component of the external potential is 
4> e (r,uj) = —eS e y, where — e is the electron charge. In the basis (§) the interband 
matrix elements of the external potential {4>e(^)}t. m = I dyi) n k{y)4>e( , r,Lj)'ip m k(y) are 
real and symmetric. Since the Hartree potential is local the matrix elements of the self- 
consistent potential are also symmetric with respect to the transverse quantum numbers 
{(ps(^)}n, m = {4>s(^)}m,n- This is also the case for local-density approximations, but 
not for the Hartree-Fock approximation, since the Fock potential is nonlocal [[D]]. The 
equation for the self-consistent potential may then be written as an eigenvalue problem 
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where e k n m = t k n — e 1 ^ are the quasi-particle excitation energies and fnm = fin ~ fn are 
the differences in occupation between initial and final states. For convenience we have 
introduced the variable {r)(u)}* fm = {^(u^^f^e^J 1 / 2 /(ti 2 u 2 - «J 2 ). The 
resonances in the FIR spectrum are the eigenvalues of the symmetric matrix 

\k;k' _ i k \2rk;k' _ C )}J^' ( f fc c k f k ' J*' \ 1/2 (q\ 

n,m;n',m' \^n,m) u n,m;n\m' n,m;n' ,m' \J n,m^n,mJ n' ,m' n' ,m' I ' \°) 

The power absorption may be found from the Joule heating 

P(cu) = - 1 dr Re {6J(f, cu)E* s (f, u)} , (9) 



where 5J(r,u) is the induced current in the wire and E s is the self-consistent electric 
field. It is 
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The Coulomb interaction may shift the poles of the matrix \j]{ui)}^ im from the Hartree 
single-particle excitations e^ m as can be seen from the eigenequation (7). The new 
poles will show up as resonances in the FIR absorption (10). Here a phenomenological 
broadening should be inserted to give the lifetime of the excitations. In a Fermi liquid 
the scattering cross-section decreases for lower energy one-particle excitations. Hence, 
one should expect that at least for low-energy one-particle excitations this broadening 
is small. The power absorption may be expanded in the eigenf unctions of the matrix 
A. The oscillator strengths can then be found and the whole FIR spectrum determined. 
Since the Hartree- RPA is a conserving approximation p2| , the longitudinal f-sum rule for 
the oscillator strengths is satisfied for arbitrary electron-electron interaction strengths. 
We study quantum wires with confining potentials of the form 
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where a and b are parameters that determine the higher order deviations from parabolic 
confinement. By the method described above, we are able to find the FIR spectra arising 
from this confining potential. 



3. Numerical results 



The functional basis of the Hartree ground state has been chosen large enough so that a 
further expansion or iteration of the Hartree-equations does not result in visual changes 
to the single-particle energy spectra or the electron density n s (f) of the ground state. To 
attain sufficient accuracy in the calculation of the absorption the size of the functional 
basis for the excited states has been chosen such that further refinement results in 
changes to the location of the absorption peaks that are smaller than a typical linewidth 
in experiments 0.1 meV. 

For the calculations we employ the usual GaAs parameters, m* = 0.067m , k = 12.4, 
where mo is the the free electron mass. The absorption is only weakly dependent on T 
in the temperature range T < 4 K for the parameters we use. The calculations have 
thus been performed for T = 1.0 K. 

The FIR spectra have been calculated in the range of magnetic fields B = — 3T. For 
a strict parabolic confinement, a = and b = 0, we have checked that the generalised 
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Kohns theorem is satisfied with a high degree of accuracy. We consider pure fourth 
order deviations, a = 0.01 to a = 0.40, and pure sixth order deviations, b = 0.001 
to b = 0.030. Mixed deviations of fourth and sixth order do not give significant new 
qualitative information. 

For small deviations, with either a / or I) / 0, a single anticrossing close to the 
2u c line appears. The position of the anticrossing as a function of the deviation a or 
b is shown in Fig. [I]. The horizontal line shows where Q p = 2u c , i.e. where a possible 
resonance between the magnetoplasmon and the cyclotron frequency should appear. The 
calculated anticrossing is to the right of the 2u; c -line for small deviations, but is shifted 
quite strongly to the left with increasing deviations. It is due to the interaction of the 
collective oscillations, the magnetoplasmons, with the first harmonic of the cyclotron 
resonance, 2u c , a so called Bernstein mode||. 

The absorption peaks as a function of the magnetic field for increasing fourth 
order deviation, a=0.03, 0.10, 0.20, 0.40, are shown in Fig. ^| and for increasing sixth 
order deviation, 6=0.003, 0.005, 0.010, 0.030, in Fig. [3[ Only absorption peaks with 
an oscillator strength more than 1% of the total oscillator strength are shown. The 
absorption peaks are labeled strong, weak, very weak and extremely weak according 
to their relative strength for a given magnetic field. Only excitation energies in the 
range 3.5-8.5 meV are shown. The location of the strong resonance at zero magnetic 
field increases in energy with increasing deviation, because of the stronger confinement. 
However for intermediate magnetic fields the strong resonance develops a slightly 
negative slope with respect to the magnetic field. This softening of the mode increases 
with increasing deviations. For large deviations the simple splitting develops into a 
complex one, see Fig. ^ for a fourth order deviation and Fig. ^ for a sixth order deviation. 
We see that strong deviations from a parabolic confinement leads to interactions with 
higher cyclotron harmonics, nu c , where n is larger than 2. The new modes may also 
interact with each other. However it is not possible from the spectra to say which 
harmonics are involved, since several harmonics can be excited in a complicated way, 
and the anticrossing occurs at or in between the nu c lines. The spectrum is now very 
rich, because of the strong coupling of the collective motion to the internal motion in 
the wire. Splittings may also be seen on the weak resonance branches for a = 0.40 in 
Fig. | and for b = 0.030 in Fig. |. 

The intensity as a function of the magnetic field and the supplied energy from the 
external potential is showed in Fig. ^ for a pure fourth order deviation, a = 0.20. Only 
resonances that have a strength of more than 1% are shown. We see that for this 
deviation there is mainly one anticrossing, all other resonances are weak. The softening 
of the modes occurs as long as two main resonances are seen, i.e as long as the modes 
are coupled. 

The effects of the deviation is larger for a higher density of the electrons, since the 
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chemical potential is higher and the 2DEG occupies states farther from the center of 
the wire that are more affected by the fourth and sixth order terms of the confinement. 

4. Summary 

The FIR spectrum for a wire with a confining potential deviating from a parabolic one 
was calculated in the time-dependent Hartree approximation. The anticrossing due to 
interaction of the magnetoplasmon with the first harmonic of the cyclotron resonance 
shifts strongly from the right of the 2u c intersection with the plasmon dispersion to the 
left with increasing deviations and the simple splitting develops into a more complex 
one, including interactions with higher harmonics of the cyclotron resonance. 

In the quantum wire huj c is not a characteristic single electron excitation energy, but 
the width of the system defines a characteristic wave vector which does not vanish with 
a decreasing deviation from the parabolic confinement. An external perturbation does 
excite plasma waves with a broad range of wave vectors, but the system responds in the 
strongest fashion for integer multiples of the 'natural wave vector'. In a homogeneous 
two-dimensional electron system the response is always with the same frequency and 
wave vector as the external perturbation, thus the asymptotic behaviour of the split- 
off modes is simple. Here the asymptotic behaviour is complicated by the fact that 
the split-off modes again interact with higher order plasmons. A detailed study of this 
phenomena will be published elsewhere. 
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Figure captions 



Figure 1. The location of the first anticrossing as a function of the magnetic field B 
and the fourth order deviation a ^ 0, b = (O) or the sixth order deviation a = 0, b =/= 
(+). The horizontal line shows where Q p = 2oj c . 30 states are occupied, L x = 240 nm, 
fiujQ = 3.94 mcV, T = 1.0 K, m* = 0.067m , n = 12.4. 



Figure 2. The location of the absorption peaks as a function of the magnetic field 
B and the energy E — Huj supplied by <p e for a pure fourth order deviation. Only 
absorption peaks with an oscillator strength larger than 1% of the total oscillator 
strength are shown. The straight lines represent the harmonics of the cyclotron 
frequency, nxuj c (n=2,3,4,5). Strong absorption is marked by O, weak one by +, very 
weak one by □, and extremely weak one by x. 30 states are occupied, L x = 240 nm, 
huj = 3.94 mcV, T = 1.0 K, m* = 0.067m , k = 12.4. 



Figure 3. The location of the absorption peaks as a function of the magnetic field 
B and the energy E = Tilu supplied by <j) e for a pure sixth order deviation. Only 
absorption peaks with an oscillator strength larger than 1 % of the total oscillator 
strength are shown. The lines are the possible interaction of the magnetoplasmon with 
harmonics of the cyclotron frequency, n x ui c (n=2,3). Strong absorption is marked 
by O, weak one by +, and very weak one by □. 30 states are occupied, L x = 240 nm, 
Tiujq = 3.94 mcV, T = 1.0 K, m* = 0.067m , k = 12 A. 



Figure 4. The intensity of the absorption peaks as a function of the magnetic field 
B and the energy E = hu> supplied by <f> e for a pure fourth order deviation a = 0.20. 
Only absorption peaks with an oscillator strength larger than 1 % of the total oscillator 
strength are shown. 30 states are occupied, L x = 240 nm, Tluiq = 3.94 meV, T = 1.0 K, 
m * = 0.067m , k = 12.4. 
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